Quantum reading under a local energy constraint 
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Nonclassical states of light play a central role in many quantum information protocols. Very 
recently, their quantum features have been exploited to improve the readout of information from 
digital memories, modeled as arrays of microscopic beam splitters [S. Pirandola, Phys. Rev. Lett. 
106, 090504 (2011)]. In this model of "quantum reading", a nonclassical source of light with 
Einstein-Podolski- Rosen correlations has been proven to retrieve more information than any classical 
source. In particular, the quantum-classical comparison has been performed under a global energy 
constraint, i.e., by fixing the mean total number of photons irradiated over each memory cell. In 
this paper we provide an alternative analysis which is based on a local energy constraint, meaning 
that we fix the mean number of photons per signal mode irradiated over the memory cell. Under 
this assumption, we investigate the critical number of signal modes after which a nonclassical source 
of light is able to beat any classical source irradiating the same number of signals. 

PACS numbers: 03.67.-a, 03.65.Ud, 42.50.-p, 89.20.Ff, 89.70.Cf 



I. INTRODUCTION 

Quantum information has disclosed a modern ap- 
proach to both quantum mechanics and information 
theory [1;]. Very recently, this field has been further 
developed into the so-called "continuous variable" do- 
main, where information is encoded and processed by 
using quantum systems with infinite-dimensional Hilbert 
spaces 0t3 (see also the recent review 6] ) . The most im- 
portant example of these systems are the bosonic modes 
of the electromagnetic field, today manipulated with 
very high precision in quantum optics labs. Thus, in 
the continuous variable framework, a wide range of re- 
sults have been successfully achieved, including proto- 
cols of quantum teleportation [t'Ji?], teleportation net- 
works [l3l - [l7| . entanglement swapping [18.- .20|, q uantum 
cryptography [2l| - l30j , quantum computation |3ll - l37| and 
cluster quantum computation |38l44ll |. 

One of the key resources in quantum information is 
quantum entanglement. In the bosonic setting, quan- 
tum entanglement is usually present under the form of 
Einstein-Podolski- Rosen (EPR) correlations [i^, where 
the quadrature operators of two separate bosonic modes 
are so correlated to beat the standard quantum limit [i^ . 
The simplest source of EPR correlations is the two-mode 
squeezed vacuum (TMSV) state. In the number-ket rep- 
resentation, this state is defined by 



|0 = (coshO-^£(tanhO"|n)Jn), 



(1) 



n=0 



where ^ is the squeezing parameter and {s,i} is an ar- 
bitrary pair of bosonic modes, that we may call "signal" 
and "idler" . In particular, ^ quantifies the signal-idler en- 



tanglement and determines the mean number of photons 
sinh^^ in each mode. Since it is entangled, the TMSV 
state cannot be prepared by applying local operations 
and classical communications (LOCCs) to a couple of 
vacua |0)g ® |0)j or to any other kind of tensor product 
state. For this reason, the TMSV state cannot be ex- 
pressed as a classical mixture of coherent states \oi) ^®\I3)- 
with a and /? arbitrary complex amplitudes. In other 
words, its P-representation (44. l45j 



m , (2) 



involves a function V which is non-positive and, there- 
fore, cannot be considered as a genuine probability dis- 
tribution. For this reason, the TMSV state is a particular 
kind of "nonclassical" state. Other kinds are single-mode 
squeezed states and Fock states. By contrast a bosonic 
state is called "classical" when its P-representation is 
positive, meaning that the state can be written as a clas- 
sical mixture of coherent states. Thus a classical source 
of light is composed by a set of m bosonic modes in a 
state 



p — I (f'ai ■ ■ ■ 



j SamV^ai,--- ,ara) ®t=i\<^k) {ak\ , 

(3) 

where V is positive and normalized to 1 . Typically, classi- 
cal sources are just made by a collection of coherent states 
with amplitudes {ai, • • • ,am}, i-e., p = 'S5'k=i l^fe) (<5fe| 
which corresponds to having 



(4) 



k=l 
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In other situations, where the sources are particularly 
chaotic, they are better described by a collection of ther- 
mal states with mean photon numbers {ni, • • • ,nm}, so 
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that 



exp(- \ak\^nk) 



TTUk 



(5) 



More generaUy, we can have classical states which are not 
just tensor products but they have (classical) correlations 
among different bosonic modes. 

The comparison between classical and nonclassical 
states has clearly triggered a lot of interest. The main 
idea is to compare the use of a candidate nonclassical 
state, like the EPR state, with all the classical states for 
specific information tasks. One of these tasks is the de- 
tection of low-reflectivity objects in far target regions un- 
der the condition of extremely low signal-to-noise ratios. 
This scenario has been called "quantum illumination" 
and has been investigated in a series of papers 46-51]. 

More recently, EPR correlations have been exploited 
for a completely different task in a completely different 
regime of parameters. In the model of "quantum read- 
ing" ^], EPR correlations have been used to retrieve 
information from digital memories which are reminiscent 
of today's optical disks, such as CDs and DVDs. A dig- 
ital memory can in fact be modelled as a sequence of 
cells corresponding to beam splitters with two possible 
reflectivities tq and ri (used to encode a bit of informa- 
tion). By fixing the mean total number of photons N 
irradiated over each memory cell, it is possible to show 
that a non-classical source of light with EPR correlations 
retrieves more information than any classical source [s^ . 
In general, the improvement is found in the regime of few 
photons {N — 1 — 100) and for memories with high reflec- 
tivities, as typical for optical memories. In this regime, 
the gain of information given by quantum reading can be 
dramatic, i.e., close to 1 bit for each bit of the memory. 
Further studies on quantum reading of memories have 
been pursued by several authors |53l459l | . In particular, 
Ref. [531 has shown that other non-classical states, such 
as Fock states, can have remarkable advantages over clas- 
sical sources. Ref. [HI] has proposed an alternative model 
of quantum reading based on a binary phase encoding. 
Ref. [55] has further studied the problem of binary dis- 
crimination in optical devices. Ref. [s^ has both pro- 
posed and experimentally implemented a model of un- 
ambiguous quantum reading. Ref. [57[ has defined the 



notion of quantum reading capacity, a quantity which has 
been also investigated in Ref. (58,]. Finally, Ref. [5^ has 
proposed explicit capacity-achieving receivers for quan- 
tum reading. 

It is fundamental to remark that an important point in 
the study of Ref. [12] is that the quantum-classical com- 
parison is performed under a global energy constraint, i.e., 
by fixing the total average number of photons N which 
are irradiated over each memory cell [see Fig. [IJa)]. Un- 
der this assumption, it is possible to construct an EPR 
transmitter, made by a suitable number of TMSV states, 
which is able to outperform any classical source com- 
posed by any number of modes. 




FIG. 1: Inset (a). Quantum reading of Ref. [S^] is formulated 
under a global energy constraint. This means that we fix 
the total average number of photons A'' irradiated over the 
memory cell. Thus, if the number of input signals is M, each 
one has an average of N/M photons, which goes to zero for 
M oo. Inset (b). In this paper, we consider an alternative 
model of quantum reading under a local energy constraint 
(locally-constrained quantum reading). In this case, we fix the 
average number of photons Ns for each input signal. Since 
the total number of signals M can be arbitrary, we have that 
the total energy irradiated over the cell MNs is generally 
unbounded. 



In this paper we consider a different kind of compar- 
ison: we fix the number of signal modes irradiated over 
the target cell (M) and the mean number of photons per 
signal mode {Ns). Under these assumptions, we com- 
pare an EPR transmitter with a classical source. Then, 
for fixed Ns, we determine the critical number of sig- 
nal modes Af after which an EPR transmitter with 
M > M^^^-' is able to beat any classical source (with 
the same number of signals M). Since we are here fix- 
ing the average number of photons per signal mode, our 
energy constraint is now local: it restricts the energy of 
each signal mode but not the energy of the total set of 
signal modes. We call this alternative model "locally- 
constrained quantum reading" [see Fig. [TJb)]. 

The difference between global and local energy con- 
straints is also discussed in Ref. G] for the general prob- 
lem of Gaussian channel discrimination. Mathematically 
speaking, both these energy constraints make the prob- 
lem of channel discrimination non-trivial in the contin- 
uous variable setting, where the use of infinite energy 
always allows one to distinguish two Gaussian channels 
in a perfect way. In the presence of a global energy con- 
straint, the error probability in the Gaussian channel dis- 
crimination is generally different from zero and the prob- 
lem is to find the minimum value. In the presence of a 
local energy constraint, the error probability goes to zero 
with the number M of signals and the general problem 
is to study its convergence, i.e., finding the best error 
exponent [6| . From this point of view, the present paper 
shows that the best convergence of the error probability 
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has to be found within the set of non-classical states. 

The remainder of the paper is structured as follows. In 
Sec. |IT] we review the basic readout mechanism of quan- 
tum reading specifying the analysis to the case of a local 
energy constraint. Then, in Sec. IIII[ we explicitly show 
how EPR correlations can be used to beat any classical 
source of light in the readout of information. Finally, 
Sec. IIVI is for conclusions. Note that we also provide 
two appendices. In Appendix [X\ we discuss the general 
mathematical methods used in our derivations, and Ap- 
pendix |B] contains some technical proofs. 



II. READOUT MECHANISM 

Here we briefly review the basic readout mechanism of 
Ref. 52] , specifying the study to the case of a local energy 
constraint. Consider a model of a digital optical memory 
(or disk) where the memory cells are beam splitter mir- 
rors with different reflectivities r — ro,ri (with ri > tq). 
In particular, the bit-value u = is encoded in a lower- 
reflectivity mirror (r = tq), that we may call a pit, while 
the bit-value u = 1 is encoded in a higher-reflectivity 
mirror (r = ri), that we may call a land (see Fig. [5]). 
Close to the disk, a reader aims to retrieve the value of 
the bit u which is stored in each memory cell. For this 
purpose, the reader exploits a transmitter (to probe a 
target cell) and a receiver (to measure the correspond- 
ing output). In general, the transmitter consists of two 
quantum systems, called signal S and idler /, respec- 
tively. The signal system S* is a set of M bosonic modes 
which are directly shined on the target cell. The mean 
total number of photons of this system is simply given 
by A^ = MNs, where Ns is the mean number of photons 
per signal mode (simply called "energy", hereinbelow). 
At the output of the cell, the reflected system R is com- 
bined with the idler system I, which is a supplementary 
set of bosonic modes whose number L can be completely 
arbitrary. Both the systems R and / are finally measured 
by the receiver (see Fig. [2]). 

We assume that Alice's apparatus is very close to the 
disk, so that no significant source of noise is present in 
the gap between the disk and the decoder. However, we 
assume that non-negligible noise comes from the thermal 
bath present at the other side of the disk. This bath 
generally describes stray photons, transmitted by previ- 
ous cells and bouncing back to hit the next ones. For 
this reason, the reflected system R combines the signal 
system S with a bath system B of M modes. These en- 
vironmental modes are assumed in a tensor product of 
thermal states, each one with Nb mean photons (white 
thermal noise). Thus, in this model we identify flve basic 
parameters: the reflectivities of the memory {ro, ri}, the 
temperature of the bath Nb, and the profile of the signal 
{M, Ns}, which is given by the number of signals M and 
the energy Ns. 

In general, for a fixed input state p at the transmitter 
(systems S, I) , Alice will get two possible output states 
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FIG. 2: Model of memory. Digital information is stored 
in a disk whose memory cells are beam splitter mirrors with 
different reflectivities: r — ro encoding bit-value u — and 
r = ri encoding bit-value u — 1. Readout. A reader is gen- 
erally composed by a transmitter and a receiver. It retrieves 
a stored bit by probing a memory cell with a signal system S 
(composed of M bosonic modes) and detecting the reflected 
system R together with an idler system / (composed of L 
bosonic modes). In general, the output system R combines 
the signal system S with a bath system B (M bosonic modes 
in thermal states). The transmitter is in a state p which 
can be classical (a classical transmitter) or non-classical (a 
quantum transmitter). In particular, we consider a quantum 
transmitter with EPR correlations between the signal and 
idler systems. In this paper, the quantum-classical compari- 
son is performed under a local energy constraint, i.e., by fixing 
the average number of photons Ns per signal mode (the signal 
system S has a total average number of photons A'^ = MNs 
which is generally unbounded). 



(To and (Ti at the receiver (systems R,I)- These output 
states are the effect of two different quantum channels, 
So and £i, which depend on the bit u = 0,1 stored in the 
target cell. In particular, we have 



CT« = {£u®1-){p) , 



(6) 



where the conditional channel acts on the signal sys- 
tem, while the identity channel I acts on the idler sys- 
tem. More precisely, we have £u — Ti-u^^ j where TZu is 
a one-mode lossy channel with conditional loss r„ and 
fixed thermal noise Nb- Now, the minimum error prob- 
ability Perr affecting the decoding of u is just the error 
probability affecting the statistical discrimination of the 
two output states, ao and cti, via an optimal receiver. 
This quantity is equal to 



[l-Diao,ai)]/2 , 



(7) 



where DjcrQ, ai) is the trace distance between ctq and 
(Ti [6Ch62| |. Clearly, the value of Perr determines the 
average amount of information which is decoded for each 
bit stored in the memory. This quantity is equal to 



where 



J = 1 - HiP,rr) , (8) 



H{x):=-x\og2X-{l-x)\og2{^-x) (9) 



is the usual formula for the binary Shannon entropy. In 
the following, we compare the performance of decoding in 
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two paradigmatic situations, one where the transmitter 
is described by a non-classical state (a quantum trans- 
mitter) and one where the transmitter is in a classical 
state (a classical transmitter). In particular, we show 
how a quantum transmitter with EPR correlations (an 
EPR transmitter) is able to outperform classical trans- 
mitters. The quantum-classical comparison is performed 
for a fixed signal profile {M, Ns}- Then, for various fixed 
values of the energy Ns (local energy constraint), we 
study the critical number of signal modes M^^^^ after 
which an EPR transmitter (with M > M*^^^^ signals) 
is able to beat any classical transmitter (with the same 
number of signals M). 



III. QUANTUM-CLASSICAL COMPARISON 

First let us consider a classical transmitter. A classical 
transmitter with M signals and L idlers is described by a 
classical state p as specified by Eq. ([3]) with m — M + L. 
In other words it is a probabilistic mixture of multi-mode 
coherent states (8)^1"^^ |a/c) (afc|. Given this transmitter, 
we consider the corresponding error probability Pg/"'''* 
which affects the readout of the memory. Remarkably, 
this error probability is lower-bounded by a quantity 
which depends on the signal profile {M,Ns}, but not 
from the number L of the idlers and the explicit expres- 
sion of the P-function. In fact, we have 52] 



pc 



>C{M,Ns) 



1 - y/l - F[Ns) 



M 



(10) 



where F{Ns) is the fidelity between 7^o(|iVy^)(iVy^|) 



and n^{\Nl'^){Nl'^\ 



the two possible outputs of 



the single-mode coherent state \Ny^){Ny^\ (see Ap- 
pendix |^ for more details). As a consequence, all the 
classical transmitters with signal profile {Af, Ns} retrieve 
an information which is upper-bounded by 



J C\ 



l-H[C{M,Ns)] 



(11) 



Now, let us construct a transmitter having the same 
signal profile {M, A^s}, but possessing EPR correlations 
between signals and idlers. This is realized by taking 
M identical copies of a TMSV state, i.e., p = |0 
where A^5 — sinh^^. Given this transmitter, we consider 
the corresponding error probability P^"^""* affecting the 
readout of the memory. This quantity is upper-bounded 
by the quantum Chernoff bound 63 67] 

Prr* < Q{M, Ns) \ [Q{Ns)f , (12) 



where 



and 



Q{Ns):^ inf Tr(0g0j-^) , (13) 
se(o,i) 



9„:=(7^„®I)(|0(el) • (14) 



Since and 9i are Gaussian states, we can write out 
their symplectic decompositions ^1] and compute the 
quantum Chernoff bound using the formula for multi- 
mode Gaussian states given in Ref. [g^] (see Appendix |X] 
for more details). Then, we can easily compute a lower 
bound 



quant 



l~H[Q[M,Ns)] 



(15) 



for the information which is decoded via this quantum 
transmitter. 

In order to show an improvement with respect to the 
classical case, it is sufficient to prove the positivity of the 
"information gain" 



G ~ J, 



quant 



J, 



class 



(16) 



This quantity is in fact a lower bound for the average 
information which is gained by using the EPR quan- 
tum transmitter over any classical transmitter. Roughly 
speaking, the value of G estimates the minimum informa- 
tion which is gained by the quantum readout for each bit 
of the memory. In general, G is a function of all the basic 
parameters of the model, i.e., G = G{M, Ns, ro,ri, Nb)- 
Numerically, we can easily find signal profiles {M,Ns}, 
classical memories {ro,ri}, and thermal baths Nb, for 
which we have the quantum effect G > 0. Some of these 
values are reported in the following table. 



M 


Ns 




ri 


Nb 


G (bits) 


1 


3.5 


0.5 


0.95 


0.01 


6.2 X 10-3 


10 


1 


0.2 


0.8 


0.01 


3.4 X 10-2 


30 


1 


0.38 


0.85 


1 


1.2 X 10-3 


100 


0.1 


0.25 


0.85 


0.01 


5.9 X 10-2 


200 


0.1 


0.6 


0.95 


0.01 


0.22 


2 X 10^ 


0.01 


0.995 


1 





0.99 



Note that we can find choices of parameters where G ~ 1, 
i.e., the classical readout of the memory does not decode 
any information whereas the quantum readout is able to 
retrieve all of it. As shown in the last row of the table, 
this situation can occur when both the reflectivities of 
the memory are very close to 1. From the first row of 
the table, we can observe another remarkable fact: for a 
land- reflectivity ri sufficiently close to 1, one signal with 
few photons can give a positive gain. In other words, the 
use of a single, but sufficiently entangled, TMSV state 
1^) (^1 can outperform any classical transmitter, which 
uses one signal mode with the same energy (and poten- 
tially infinite idler modes). 

Here an importat point to remark is that, once that we 
find a positive gain G > 0, this positivity is preserved if 
we increase the number of signals M. In other words, if G 
is positive for some M, then it is positive for every M > 
M (keeping the other parameters fixed.) This is trivial 
to prove. In fact, G{M) > is equivalent to Q{M, Ns) < 
C{M, Ns) which is equivalent to 



(17) 
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according to Eqs. pO)) and (IT2]) . This means that 

Q < (^1 - V 1 - F*^ ^ . (18) 

For every M > M, we then have 

Q^' < (1 - VT^)" , (19) 

where m := M/M and x :— F*^. Now, we can use the 
algebraic inequahty 

(1 - \/]~^)™ < 1 - VI -x™ , (20) 
which holds for every m > 1 and x Cz [0, 1]. Then, we get 
qM < 1 _ ^ (21) 

which is equivalent to 

Q{M,Ns)<C{M,Ns) , (22) 

for every M > M. 

Thanks to this property, for given reflectivities {ro, ri} 
and bath temperature iV^, i.e., for a fixed memory, if a 
quantum transmitter with signal profile {M, Ns} outper- 
forms the classical transmitters, then any other quantum 
transmitter with the same energy Ng and M > M is also 
able to beat the classical readout. 



A. Ideal memories 

According to our numerical investigation, the quan- 
tum readout is generally more powerful when the land- 
reflectivity is sufficiently high (i.e., ri > 0.8). For this 
reason, it is very important to analyze the scenario in the 
limit of ideal land-reflectivity (ri = 1). Using the termi- 
nology of Ref. [13 , we call an "ideal memory" a classical 
memory with ri = 1. Clearly, this memory is completely 
characterized by the value of its pit-reflectivity rg. For 
ideal memories, the quantum Chernoff bound of Eq. ([12]) 
takes an analytical form given by the "Chernoff term" 



^' [l + {l-^)Ns]'' + NB{2Ns + l){l-ro)' 

(23) 

and the classical bound of Eq. ([TU)) can be computed 
using 

FiNs) = exp[-7-i(l - ^^Ns] , (24) 

where 7 := 1 -I- (1 — ro)NB (see Appendix for more 
details). Using these formulas, we can study the behav- 
ior of the gain G in terms of the remaining parameters 
{M, Ns,ro, Nb}- Let us consider an ideal memory with 
generic tq € [0,1) in a generic thermal bath Nb > 0. 
For a fixed energy Ns, we consider the minimum num- 
ber of signals M^^®^ above which G > 0. To be precise, 



the critical number M^^'=^^ that we consider is a solution 
of the equation G = 0. From this real value we derive 
the minimum number of signals (which is an integer) by 
taking its ceiling function [M*^^^)]. The critical number 
can be defined independently from the thermal 
noise Nb by performing a numerical maximization over 
Nb- Then, for a given value of the energy Ns, the crit- 
ical number M^-^^^ becomes a function of tq alone, i.e., 
^/(Ws) = M^^s\ro). Its behavior is shown in Fig. |3] for 
different values of the energy. 




0.2 0.4 0.6 0.8 1 

FIG. 3; Number of signals M (logarithmic scale) versus pit- 
reflectivity ro. The curves refer to Ns = 0.01, 0.1 and 0.5 
photons. For each value of the energy Ns, we plot the critical 
number M^^^'(ro) as function of ro. All the curves have an 
asymptote at ro = 1. For A''s > 2.5 photons (curves not 
shown), we have another asymptote at ro = 0. 

It is remarkable that, for low-energy signals (Ns — 
0.01 — 1 photons), the critical number M'^^^)(ro) is fi- 
nite for every tq € [0,1). This means that, for ideal 
memories and low-energy signals, there always exists a 
finite number of signals M^^^-* above which the quantum 
readout of the memory is more efficient than its classical 
readout. In other words, there is an EPR transmitter 
with M > M (^^^ able to beat any classical transmit- 
ter with the same number of signals M. In the low- 
energy regime considered, M'-^^^(ro) is relatively small 
for almost all the values of tq, except for rg — > 1 where 
Af(^^^(ro) 00. In fact, for tq ~ 1, we derive 

M(^^)(ro) ~ [msi2Ns + 1)(1 - roT' , (25) 

which diverges at tq = 1. Such a divergence is expected, 
since we must have Pg'^"""* — P^l'^^'^ = 1/2 for ro = ri 
(see Appendix [B] for details). Apart from the divergence 
at To ~ 1, in all the other points tq G [0, 1), the crit- 
ical number M'-^^^(ro) decreases for increasing energy 
Ns (see Fig. [3]). In particular, for Ns — 1 photon, we 
have M(^®'(ro) ~ 1 for most of the reflectivities tq. In 
other words, for energies around one photon, a single 
TMSV state is sufficient to provide a positive gain for 
most of the ideal memories. However, the decreasing 
trend of M(^^)(ro) does not continue for higher energies 
{Ns > 1). In fact, just after TVs = 1, M'^^^\ro) starts 
to increase around ro = 0. In particular, for Ns > 1, we 
can derive 

M(^^)(0) ~ (ln2)[21n(l + A^s) -iVs]-^ , (26) 
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which is increasing in Ns, and becomes infinite at Ns — 
2.5. As a consequence, for Ng > 2.5 photons, we have a 
second asymptote appearing at tq = (see Appendix IbI 
for more details). This means that the use of high-energy 
signals {Ns > 2.5) does not assure positive gains for 
memories with extremal reflectivities ro — and ri — 1. 



IV. CONCLUSION 

In conclusion, we have considered the basic model of 
digital memory studied in Ref. [s^] , which is composed of 
beam splitter mirrors with different reflectivities. Adopt- 
ing this model, we have compared an EPR transmit- 
ter with classical sources for fixed signal profiles, finding 
positive information gains for memories with high land- 
reflectivities (ri > 0.8). Analytical results can be derived 
in the limit of ideal land-reflectivity (ri = 1) which de- 
fines the regime of ideal memories. In this case, by fixing 
the mean number of photons per signal mode (local en- 
ergy constraint), we have computed the critical number 
of signals above which an EPR transmitter gives positive 
information gains, therefore beating any classical trans- 
mitter. For low-energy signals (0.01 — 1 photons) this 
critical number is finite and relatively small for every 
ideal memory. In particular, an EPR transmitter with 
one TMSV state can be sufficient to achieve positive in- 
formation gains for almost all the ideal memories. 

Thus our results corroborate the outcomes of Ref. [l^l 
providing an alternative study which considers a local 
energy constraint instead of a global one. As discussed 
in Ref. 's^, potential applications are in the technology 
of optical digital memories where we could increase data- 
transfer rates and storage capacities. Another potential 
application is the readout of organic digital memories 
which are extremely photosensitive at high frequencies. 
In this case, the use of faint narrowband quantum signals 
could safely read the data without damaging the storing 
devices. 



Appendix A: Methods 

Here we provide more details about the methods used 
in our computations. We start with a brief review of the 
bosonic Gaussian states (Appendix lA 1|) . Then, we give 
a detailed description of the readout problem discussing 
the main techniques for computing the performances of 
classical and nonclassical transmitters (Appendix lA 2| . 
Finally, we consider the special case of ideal memories, 
for which we can derive simple analytical formulas (Ap- 
pendix IA3p . 



1. Bosonic systems and Gaussian states 

A bosonic system is generally composed of n modes. 
This means that the system is associated to a tensor 



product Hilbert space 'H®" and described by a vector 
of quadrature operators 



which satisfy the commutation relations [xk, xi] 
where 



e 

1=1 



1 

-1 



(Al) 
(A2) 



defines a symplectic form (correspondingly, a real matrix 
S is called "symplectic" if SHS^ = fl). 

By definition, a bosonic state p is "Gaussian" if its 
Wigner function is Gaussian f|. As a result, a Gaussian 
state p is fully characterized by its first and second order 
statistical moments. These are the displacement vector 

X := Tr(xp) , (A3) 

and the covariancc matrix (CM) V, with generic element 

Vki := ^Tr {{xk, xijp) - XkXi , (A4) 

where {, } is the anticommutator. The CM is a 2ri x 2n 
real symmetric matrix which must satisfy the uncertainty 
principle [gI. |69| 



V + in>o 



(A5) 



According to Williamson's theorem every CM V can 
be decomposed in the form 



V = SWS^ , 

where S is a symplectic matrix and 

/-I 



(A6) 



W 



i=l 



\ 



(A7) 



is called the "Williamson form" of V. In this matrix, the 
diagonal elements \v\.--- , i^„} represent the "symplec- 
tic spectrum" of V. The symplectic spectrum provides 
powerful ways to express physical properties of the Gaus- 
sian state. For instance, the uncertainty principle can be 
formulated as [1, [6l| 



V > , Vi>\ 



(A8) 



2. Quantum reading versus classical reading 

In our model of classical memory, each memory cell 
is represented by a beam splitter mirror with two possi- 
ble reflectivities, i.e., the pit- reflectivity tq and the land- 
reflectivity r\. This dichotomic choice r € {tq, ri} is used 
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to encode an information bit u € {0, 1} in the memory 
cell. Then, one side of the memory is subject to decoding, 
while the other side is affected by white thermal noise, 
with average photon number per mode equal to Nb- It is 
clear that this model of memory corresponds to a prob- 
lem of Gaussian channel discrimination. In fact, each 
memory cell can be seen as an attenuator channel, trans- 
forming an input signal mode into an output reflected 
mode. In particular, this attenuator channel has a trans- 
mission efficiency (or "linear loss" ) which is given by the 
dichotomic reflectivity of the cell r g {vq, ri}, and a ther- 
mal noise which is fixed and equal to Nb- Depending on 
the bit u G {0, 1} which is stored in the cell, we then have 
two possible attenuator channels, that we denote by TZo 
and TZi. In other words, the unknown bit u stored in the 
cell is encoded into a conditional attenuator channel 7?.„. 

Let us describe explicitly the action of TZu- Given the 
quadratures := {qs,Ps) of an input signal mode s, the 
quadratures x^ of the output reflected mode r are given 
by the Heisenberg relation 



= ^/^Xs + Vl - 



(A9) 



where r € {?'Oi''i} and Xf, are the quadratures of a bath 
mode b. In particular, the bath mode is described by a 
thermal state pbiNs) with Nb average photons, i.e., a 
Gaussian state with zero mean and CM 



Vfc = {2Nb + 1)1 



(AlO) 



Once that we have specified the action of a memory cell 
over an arbitrary signal mode s, we can analyze its full 
action on an arbitrary transmitter. In general, we have 
a system S of M signal modes impinging on the cell, 
besides an ancillary system I of L idler modes which 
bypass the cell. At the output of the cell, the system 
R of the M reflected modes is combined with the idler 
system / in a joint measurement at the receiver. The 
fundamental parameters of the transmitter are contained 
in its signal profile {M,Ns}, which is composed by the 
number of signal modes AI and the average number of 
photons per signal iVg. 

Let us denote by psi the global state of the input sys- 
tems {S,I}. The memory cell does not affect the idler 
system J, but acts on the signal system S by coupling ev- 
ery signal mode s G S with an independent thermal mode 
6, which belongs to a bath system B in the multiniode 
thermal state 



Pb = Pb{NB) 



Since the action on the signal system S is one-mode and 
conditional, the global state of the output systems {R, 1} 
can be written as 



PRiiu) = (7^: 



{p. 



SI) 



(A12) 



where X'^^ is the identity channel acting on the idler sys- 
tem. For a fixed state psi at the transmitter, we have 
a conditional output state pri(u) at the receiver, which 



depends on the bit u stored in the memory cell. Thus, 
the minimum error probability in decoding the stored bit 
is just the error probability affecting the optimal discrim- 
ination of the two output states pr,i{0) and Pfl,/(f). As 
we know, this error probability is equal to jGOiJ 



''err = l{l-D[pRl{0),pRl{l)]} , 



(AI3) 



where Z)[/9_r/(0), p_R/(l)] is the trace distance between 
Pri{0) and piii{l). Clearly, the value of Perr determines 
the average amount of information which is decoded for 
each bit stored in the memory. This average information 
is equal to J = 1 — H{Perr), where H{x) is the usual 
formula of the binary Shannon entropy. 

In our work, we estimate the average decoded infor- 
mation J in two paradigmic situations, i.e., for a quan- 
tum transmitter with EPR correlations (Jq), and for a 
generic classical transmitter (Jc)- By fixing the signal 
profile {M,Ns}, we compare Jq and Jc- More exactly, 
we fix all the basic parameters of the model, i.e., besides 
fixing the signal profile {M, Ns}, we also fix the reflectiv- 
ities of the memory {ro,ri} and the thermal noise Nb- 
Then, we investigate what are the values of the basic 
parameters {M,Ns,ro,ri,NB} for which Jq > Jq- In 
particular, for proving this enhancement, we compare a 
lower bound of Jq with an upper bound of Jc- 



a. Classical transmitters 

Let us start considering an arbitrary classical trans- 
mitter. For a classical transmitter with M signals and L 
idlers we can exploit the classical discrimination bound 
proven in Ref. [52]. The minimum error probability 
Pg/j!''^ affecting the readout of the memory cell is lower- 
bounded by C{M,Ns) in Eq. ^ where F{Ns) is the 
fidelity between the two states 'R-oi\Nl^^){Ng^^\) and 
7^l(|7Vy^)(iVy^|). Using thejormula of the fidelity for 



single- mode Gaussian states |72h74| |. we get 



FiNs) 



where 



1 



exp 



Ns 



7 = 1 + (2 - ro - ri)iV£ 



(A14) 



(A15) 



(All) and 



e ^ ml II (l - r,)[l + (1 - r,)NB] - (Af6) 



Notice that the lower-bound C{M, Ns) depends on the 
signal profile {M,Ns}, but not from the number L of 
idlers and the explicit P-representation describing the 
classical state of the transmitter. As a consequence, all 
the classical transmitters with the same signal profile 
{M,Ns} are lower-bounded by C{M,Ns)- The average 
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information Jc which is decoded from the memory ceh 
is upper-bounded by the quantity 



J das 



l-H[C{M,Ns)] 



(A17) 



b. Quantum transmitter 



Now, let us consider a quantum transmitter with the 
same signal profile {M, iVs} but possessing EPR corre- 
lations between signals and idlers (EPR quantum trans- 
mitter) . In this case, we have the same number of signals 
and idlers [M = L), and the global state for the input 
systems {S, 1} is a tensor product of Af identical two- 
mode squeezed vacuum states, i.e.. 



(A18) 



where the single-copy state |^)^^ (^| refers to a single pair 
of signal and idler modes {s,i} G {S,I}. Recall that a 
two- mode squeezed vacuum state 1^)^,^ (^| is a Gaussian 
state with zero mean and CM 



{2Ns + 1)1 2v/iVs(iVs + l)Z \ 
^ 2^Ns{Ns + l)7. (2iVs + l)I / ' 

(A19) 

where I = diag(l,l), Z — diag(l,— 1) and the squeez- 
ing parameter ^ is connected to the signal-energy by the 
relation Ns — sinh^ ^. At the output of the cell, the 
conditional state of the systems {R, /} is given by 



where 



PRi{u)^p„{ur^'' , (A20) 



P„(u) = (7^„®T)(|0,,(^|) (A21) 



is the single-copy output state, i.e., describing a single 
pair of reflected and idler modes {r, i} G {i?, /}. In fact, 
since the memory cell corresponds to a one-mode chan- 
nel and the state of the transmitter to a tensor product, 
the output state at the receiver is also a tensor product 
state. In particular, it corresponds to M identical copies 
of the two-mode state of Eq. (|A2ip . Then, the decoding 
of u corresponds to the M-copy discrimination between 
the two states Pri(O) and Pri(l)- The corresponding min- 
imum error probability Pg'^°"* can be upper-bounded by 
the quantum Chernoff bound, i.e.. 



Prr* < Q{M,Ns) \ [Q{Ns)] 



M 



where 



Q{Ns) := inf Tv [p„{Or p„{iy-'] 

se{o,i) 



(A22) 



(A23) 



Notice that the single-copy state Pri{u) is a Gaussian 
state with zero mean and CM 



Vrt('u) 



[ruP + (1 



Z 



1)Z 



1 



(A24) 



where 

p:=2Ns + l , /3 2iVB + 1 . (A25) 

Since Pri(O) and Pri(l) are two-mode Gaussian states, we 
can compute the quantum Chernoff bound Q(M, Ns) by 
using the formula of Ref . [63] , which exploits the symplec- 
tic decomposition of the Gaussian states. It is important 
to note that the CM of Eq. (|A24p is in a special form, 
for which we can easily provide analytical expressions 
for both the symplectic spectrum and the diagonalizing 
symplectic matrix S. In fact, let us set 



a := ruP + (1 - r^)/? , b := p , 



(A26) 



and 



c := v/r„(/i2 - 1) > , (A27) 
so that the CM takes the special normal form 



V = 



al cZ 
cZ 61 



(A28) 



The corresponding symplectic eigenvalues are given by j 



and 



(A29) 



(A30) 



where y :— {a + bY — 4c^ > 4. Thus, the Williamson 
form of V is given by 



2 I i^ + b^a)I 



(A31) 



The symplectic matrix S which realizes the symplectic 
decomposition V = SWS"'" is given by the formula 



S = 




where 



x± 



la + b±y/y 



2^/y 



> 



(A32) 



(A33) 



Now, by expliciting the symplectic diagonalization for the 
two possible cases u — and u = 1, we have 

V„{u) = S{u) [i^i{u)I®i^2{u)I] S{uf . (A34) 

Using this decomposition, we can compute the quantum 
Chernoff bound by means of the formula of Ref. [67| . 
Unfortunately, the analytical expression is cumbersome, 
but we can easily derive numerical values for every choice 
of the parameters. 

It is clear that, from the upper bound < 
Q{M,Ns), we can derive a lower bound for the aver- 
age information Jq which is decoded via this quantum 
transmitter. This lower bound is simply given by 



quant 



1 - H[QiM, Ns)] . 



(A35) 
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c. Comparison 

In order to compare quantum and classical reading, we 
fix the basic parameters of the model {M, Ns, rQ,ri, Ng} 
and we consider the difference 

- Jclass , (A36) 

that we call "information gain" . It is trivial to check that 
G<Jq-Jc. (A37) 

In other words, G is a lowerbound for the average in- 
formation which is gained by using the EPR quantum 
transmitter instead of any classical transmitter. A pos- 
itive gain (G > 0) is a sufficient condition for the supe- 
riority of the quantum reading ( Jq > Jc)- In general, 
this quantity is a function of all the basic parameters of 
the model, i.e., G = G{M, Ns,ro,ri,NB)- Numerically, 
we can find signal profiles {M,Ns}, classical memories 
{ro,ri}, and thermal baths Nb, for which we have the 
quantum effect G > 0. Some of these values are shown 
by the table in the main text. 



with CM given in Eq. (jA19[) . Then, we can apply the for- 
mula of Ref. ^TlJ for the quantum fidelity between mul- 
timode Gaussian states. We achieve 

Q{Ns)^F[p„{0)AOs^ = 

1 

~ [1 + (1 - V^)Ns]^ + Nb{2Ns + 1)(1 - ro) ' 

(A41) 

which is the result given in the main text. 

In the case of a classical transmitter, we just have to 
consider the lower bound of Eq. (1101) where now we set 
ri = 1 in the expression of the fidelity given in Eq. (jA14|) . 
One can easily check that the resulting fidelity takes the 
analytical form given in Eq. (|24l) . 

Appendix B: Technical proofs 

Here, we explicitly prove the asymptotic expansions 
which have been presented in the main text and used for 
the analysis of the ideal memories. 



Ideal memories 



1. General asymptote (ro — 1) 



Quantum reading is generally more powerful when the 
land-reflectivity is sufficiently high (i.e., ri > 0.8). For 
this reason, it is important to analyze the scenario in 
the limit of ideal land- reflectivity (ri = 1), deflning the 
so-called "ideal memories". In the presence of an ideal 
memory, one of the two possible outputs of the cell is just 
the input state, i.e., we have 



/Ofl/(l) = Psi ■ 



(A38) 



Clearly, this fact leads to a simpliflcation of the calculus. 
In the case of an EPR quantum transmitter, the input 
state is pure and given by Eq. (jAlSp . As a consequence, 
we have 



(giM 



10.. (ei 



0M 



(A39) 



i.e., one of the two output states is pure. As a con- 
sequence, the quantum Chernoff bound can be reduced 
to the computation of the quantum fidelity. In fact, we 
have (nl 



QiNs) 



inf Tr [p„(0)Vr»(l)'" 
se(o,i) 



inf Tr 

se(04) 



Pr^{0y\Os^{^\ 



1-s 



lim Tr 

s-i-l- 



Pr^iOy\Os^ i^l 



(A40) 



where the fidelity F[pri{0), lO^j between a mixed 

two- mode Gaussian state pri(0) with CM given in 
Eq. (jA24|) and a pure two-mode Gaussian state 1^)^,^ (^| 



According to Fig. [3l the critical number M^^^\rQ) 
diverges for rg — ?> 1. Let us analyze the behavior of G 
around the singular point ro = 1, by setting ro = 1 — e 
and expanding G for e — > 0+. It is easy to check that, 
for every Nb , we have G > if and only if 



M > [4Ns{2Ns + l)e]"^ 



(Bl) 



In particular, in the absence of thermal noise {Nb — 0), 
we have 



G = 



MNsi-iMNs 



;ln2 



(B2) 



which is positive if and only if AI > {ANs)~^- 

These conditions are easy to prove. In fact, note that 
G > if and only if 



QiM,Ns)-C{M,Ns) <0 



(B3) 



Thus, let us expand A = A(M, Ns, iVs, 1 - e) at the first 
order in s. For a given Nb > 0, we have 



2 



M{Ne 



{MNBef^ . 

which is negative if and only if 

Nb 



Ns 



2NBNs)e 



M > 



{Nb +Ns + 2NBNsfe 
Notice that k(Nb) is maximum for 

N*s = Nsil + 2Nsr' 



<Nb) 



(B4) 



(B5) 



(B6) 
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Then, for every Nb > 0, we have A < if and only if 

1 



ANsi2Ns + l)e 



Now, let us consider the particular case of Nb 
this case, we have the first-order expansion 



(B7) 
0. In 



A 



{MNs f^^ [1 



2(MiVs)^/']e/4 + 0(e2) , (B8) 



or, equivalently, the second-order expansion of G given in 
Eq. (|B2]). It is clear that A < 0, i.e., G > 0, when M > 
1/ANs- However, this condition is less restrictive than 
the one in Eq. (jB7|) which, therefore, can be extended to 
every Nb > 0. 



2. High-energy asymptote (t-q — 0) 

Let us analyze the behavior of Af'^^^^(ro) for Ng > 1 
and ro — 0. One can check that, for Ns > 1, the greatest 
value of M(^s)(0) occurs when Nb = 0. In this case, i.e., 
for ro ~ Nb ~ and ri = 1, we have 



Q.(M,Ns) 



(1 + Ns) 



-2M 



(B9) 



and 



C{M,Ns) 



1 e 



'MNs 



(BIO) 



Let us consider the critical value M(^s)(0) of M such 
that G{M, Ns) = 0, which is equivalent to Q = C. We 
also consider the value M such that Q — C°°. We find 
that M(^^'(0) ~ M with very good approximation when 
> 1 (see Fig. |4]). Then, for every Ns > 1, we can set 



M(^s)(o) ~ M = (In 2) [21n(l + Ns) ~ Ns 



(Bll) 



The latter quantity becomes infinite for 21n(l -|- Ns) 
Ns, i.e., for TVs > 2.51 photons. 



M 




FIG. 4: Minimum number of signals M versus energy Ns- 
The solid curve represents M''^s)(o) while the dashed curve 
represents M. Notice that the minimum number of signals is 
actually given by \M] where [■ • • ] is the ceiling function. 
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